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Abstract. 

Satellite lightning imagers such as the NASA Tropical Rainfall Measuring Mission Lightning Imaging 
Sensor (TRMM/L1S) and the future GOES-R Geostationary Lightning Mapper (GLM) are designed to 
detect total lightning (ground flashes + cloud flashes). However, there is a desire to discriminate ground 
flashes from cloud flashes from the vantage point of space since this would enhance the overall 
information content of the satellite lightning data and likely improve its operational and scientific 
applications (e.g., in severe weather warning, lightning nitrogen oxides studies, and global electric circuit 
analyses). A Bayesian inversion method was previously introduced for retrieving the fraction of ground 
flashes in a set of flashes observed from a satellite lightning imager. The method employed a constrained 
mixed exponential distribution model to describe the lightning optical measurements. To obtain the 
optimum model parameters (one of which is the ground flash fraction), a scalar function was minimized 
by a numerical method. In order to improve this optimization, a Grobner basis solution was introduced to 
obtain analytic representations of the model parameters that serve as a refined initialization scheme to the 
numerical optimization. In this study, we test the efficacy of the Grobner basis initialization using actual 
lightning imager measurements and ground flash truth derived from the national lightning network. 
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1. INTRODUCTION 


3. TEST OF THE GROBNER INITIALIZATION 


When viewing lightning from space at optical wavelengths, the cloud multiple scattering medium obscures the view thereby preventing one 
from easily determining what flashes strike the ground. However, recent studies have made some progress in estimating the ground flash 
fraction , cjc in a set of N flashes observed from space [Koshak (2010), Koshak and Solakiewicz (2011), and Koshak (2011)]. Knowledge of 
dss> important for better understanding: 


• Severe Weather, 

• Lightning Nitrogen Oxides Chemistry/Climate Studies, and 

• Global Electric Circuit. 

In the study by Koshak (2011), a Bayesian inversion method 
was introduced for retrieving the fraction of ground flashes in 
a set of flashes observed from a (low earth orbiting or geostationary) 
satellite lightning imager. This method has formed the basis of a 
Ground Flash Fraction Retrieval Algorithm (GoFFRA) that is being 
tested as part of GOES-R Geostationary Lightning Mapper (GLM) 
risk reduction. 


BASIC 

QUESTION: 


2. THE BAYESIAN RETRIEVAL 

Figure 2 highlights the mathematical attributes of the Bayesian 
retrieval process. The // ' s are the population mean Maximum 
Group Areas (MGAs) decremented by the nadir footprint area of 
the space-based lightning imager employed. In summary, the 
distribution of y (nadir footprint decremented MGA) is assumed 
to follow a mixed exponential distribution; one finds the 3 para- 
meters of this distribution that maximizes a Bayesian scalar function. 



Figure 1. Example of lightning optical emissions as seen from space (STS-52 video). 



To test the Grobner Initialization approach and compare it with the Original Initialization, we perform simulated retrievals. To simulate the 
observation vector y (i.e., the set of MGAs each decremented by the nadir footprint area), we generate the components of y from a mixed 
exponential distribution having parameters (^C£ rRLJE , fi gTRLJE 9 fi c TR LJE ). These true values are compared with the Bayesian 
retrieved values to assess retrieval errors. The true values are varied over a wide range, and a total of 100 retrievals are performed in each 
ground flash fraction bin (0-0.05, 0.05-0.1, ... , 0.95-1.0) . The mean retrieval error (and std. deviation) are plotted below. 
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4. GRID SEARCH 


The Grobner Initialization + Minimization does not improve upon the Original Initialization + Minimization . Moreover, we have found that a 
simple grid search over the 3 parameters (for finding the maximum in S) can beat out both of these methods as shown in the (below left) 
retrieval error plot. Each retrieval takes only tens of seconds on a Dell T3500 Linux machine (depending on grid size) and no normal prior 
constraints were needed. Plots of the S scalar can also be generated to provide more insight (see sequence of plots below right). 


Plot of S ( /A, , ju c )for sample Grid Search retrievals; cc TRUE varies from panel to panel: 
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